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Abstract. This paper presents a non-classic approach to narrow field astrometry that offers a significant improvement over 
conventional techniques due to enhanced reduction of atmospheric image motion. The method is based on two key elements: 
apodization of the entrance pupil and the enhanced virtual symmetry of reference stars. Symmetrization is implemented by 
setting special weights to each reference star. Thus a reference field itself forms a virtual net filter that effectively attenuates the 
image motion spectrum. Atmospheric positional error was found to follow a power dependency A ~ R*A2£)-r/2+i/3 on an g U j ar 
field size R and aperture D\ here k is some optional even integer 2 < k < s/8N + 1 - 1 limited by a number N of reference 
stars, and p < 1 is a term dependent on k and the magnitude and sky star distribution in the field. As compared to conventional 
techniques for which k = 2, the improvement in accuracy increases by some orders. Limitations to astrometric performance of 
monopupil large ground-based telescopes are estimated. The total atmospheric and photon noise for at a 10 m telescope at good 
0.4" seeing was found to be, depending on sky star density, 10 to 60 pas per 10 min exposure in R band. For a 100 m telescope 
and FWHM=0.1" (low-order adaptive optics corrections) the potential accuracy is 0.2 to 2 pas. 
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1. Introduction 

Atmospheric turbulence affects astronomical observations with 
respect to both image quality degradation and random im¬ 
age motion. With improvement of image detectors the second 
source of errors has become dominant in astrometric measure¬ 
ments. The very slow progress in overcoming atmospheric er¬ 
rors finally turned to development of high-accurate astrome¬ 
try from space missions. The first successful Hipparcos mis¬ 
sion yielded data with an accuracy of 1 mas, the GAIA project 
(Perryman et all. 2001) is aimed at 10 pas accuracy level. 

The most promising direction of optical ground-based as¬ 
trometry is now related to the infrared long-baseline interfer¬ 
ometers; the atmospheric limit for this kind of instruments is 
about 10 pas (Shao & Colavita 1992; Paresce et all. 2002). 
Astrometry with monopupil long-focus telescopes is limited 
by much larger A ~ 1 mas atmospheric error, which was con¬ 
firmed by numerous experimental and theoretical data. Thus, 
Gatewood (1987), on the basis of measurements with the 
Multichannel Astrometric Photometer (D = 0.76 m), derived 
A =: 3 mas/hr with a reference frame of 10 - 20'. Using the 
same data, Han (1989) suggested a model of a image motion 
power spectrum according to which the distance between a 1' 
separated double star can be measured with an accuracy of 
1 mas/lhr. The best precision of 150 pas/hr ever obtained 
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with ground-based telescopes was demonstrated by Pravdo & 
Shaklan (1996) with the 5-m Palomar telescope in a field of 90 
arcsec. 

A theoretical aspect of the image motion problem was stud¬ 
ied by Lindegren (1980) who derived analytical expressions for 
relative image motion of a double star. For very narrow fields, 
A 2 ~ (fih) 1 !) 4,1 7' 1 demonstrated a rather weak power depen¬ 
dency of the variance of image motion A 2 on the double star 
angular separation p, objective diameter D and exposure T (h 
is the turbulent layer height). The above dependence shows that 
an improvement of A from 1 mas to say 10 pas will require un¬ 
realistically large D and T, unless p is limited to uselessly small 
angles. The situation is somewhat better when a very dense star 
field is used as a reference. For circular star samples of a ra¬ 
dius p the dependency becomes A 2 ~ {phf'^D 2 T~ 1 which is 
a factor (ph/D ) 2 ^ 3 better than that valid for a single reference 
star. Lazorenko (2002a), making allowance for the second or¬ 
der term of the random phase fluctuations over the pupil, has 
proved that the dependency is A 2 ~ (p/?) ll : /J 3 T ’, or a fac¬ 
tor ph/D still stronger. He also suggested a symmetrizing pro¬ 
cedure by means of which any arbitrary distributed sample of 
stars can be utilized as a virtually symmetric reference frame. 

Nevertheless even the use of symmetric (or symmetrized) 
reference fields not can provide accuracies like those expected 
from interferometric and/or space mission techniques. Thus, 
Lazorenko (2002a), considering a 10 m telescope, 1' diame- 





2 


Lazorenko P. & Lazorenko G.: Filtration of atmospheric noise 


ter reference frame and all the turbulence concentrated at 2.8 
km height, found that atmospheric noise is 5 /vas/hr at moderate 
seeing. Scaling this data to more realistic effective heights of 14 
-17 km yields A = 120 - 300 /ras/hr, which is the limit of this 
technique. An efficient solution of the problem suggested in 
the current paper is based on enhanced symmetrization of ref¬ 
erence fields, which is an improved version of the former pro¬ 
cedure (Lazorenko 2002a). This new approach allows a much 
better filtration of random wave-front distortions, reducing at¬ 
mospheric noise to 1-10 /(as/hr with a 10 m aperture. 

This paper presents the analysis of atmospheric limits to 
astrometric performance of one-aperture telescopes. We em¬ 
phasise the use of very large 10 - 100 m telescopes since we 
want to obtain limiting estimates of ground-based astrometry. 
Astrophysical drives and possible applications for microarc- 
second astrometry are not discussed here although they were 
considered in some papers (Perryman et al. 2001; Pravdo & 
Shaklan 1996; Paresce et al. 2002). In Section 2 we describe a 
model of a differential image motion spectrum. The image mo¬ 
tion spectrum filtration described in Section 3 is based on two 
key elements: apodization of the entrance pupil (Section 4) and 
the enhanced virtual symmetry of reference groups (Section 5). 
In Section 6 we consider asymptotic properties of high-order 
symmetry dense reference fields and in Section 7 we discuss 
modifications of plate reductions necessary to implement the 
new technique. Estimates of the image motion error integrated 
over the atmosphere are given in Section 8. Astrometric per¬ 
formance of very large telescopes under assumption of a two- 
component error budget is analyzed in Section 9. 


2. A model of differential image motion power 
spectrum 

Atmospheric random motion of the star image centroid is 
caused by turbulence which is known to be concentrated in a 
limited number of thin atmospheric layers. Each layer is de¬ 
scribed by its height h above the ground, a thickness A/z which 
is typically 10 - 100 m, rarely up to 800 m (e.g. Barletti et al. 
1977; Redfern 1991; Vinnichenko et al. 1976) and a refractive- 
index structure constant C 2 . Turbulent motion in each layer is 
isotropic in horizontal directions, and refractive-index fluctu¬ 
ations are described by the structure function D„( r) = Cf t r p , 
where r is a vector connecting two points in a 2-D layer’s plane, 
and p is a constant. Some experimental and theoretical studies 
discussed by Lazorenko (2002a) testify that p is to be consid¬ 
ered as a variable term with a typical range of variations 1/2 — 
2/3 which at some particular atmospheric conditions is widened 
from -1/3 to +1. In the present study, for the sake of universal¬ 
ity of the derived equations, this quantity is considered as a 
variable, but numerical estimates are made with p — 2/3 only. 

Due to the small thickness of layers A/z <s /z, the each one 
is considered as a thin horizontal phase screen that produces 
random phase distortions of the light wave front. A plane wave 
that crosses the screen thus gains some random phase <p(x, y) 
where x, y are Cartesian coordinates in the phase screen plane. 
Given the assumption of homogeneity and isotropy of turbu¬ 


lence (at least in a horizontal plane), the power spectral density 
of a phase is a function 

F*iq) = c^q-”- 3 ( 1 ) 

of a circular spatial frequency q. Here c (!) is a factor accounting 
for turbulence strength. In case of the Kolmogorov turbulence, 
cq is related to C 2 via the expression (Tatarsky 1961) 

c 0 = 0.033(27r/d) 2 (27r)“ 2/3 C;A/z (2) 


where A is wavelength. 

Perturbations of a light phase (f>(x,y ) lead to random dis¬ 
placements / of star image centroids in the focal plane. The 
magnitude of / measured in some arbitrary direction r' is pro¬ 
portional to the gradient of the phase d<p{x, y) / dr’ averaged over 
the screen section involved in image formation. For the instant 
exposure T = 0 it is a circular area of a diameter D equal to the 
telescope diameter and centered at the point r where the phase 
screen is crossed by light beams passing from the star to the 
centre of the telescope pupil. A quantity / as a function of r 
vector is given by the convolution (Martin 1987; Conan et al. 
1995) 


m = 


A dc/>(r) P( r) 
2n dr' cr 


(3) 


where P( r) is the entrance pupil function that determines opti¬ 
cal transmission of the objective and 


cr — 



P( r) dr 


(4) 


is the effective area of the entrance pupil. For a classical fully 
transparent monopupil 


P( r) = 


1, r<D/2 
0, r>D/2 


(5) 


At finite exposure T > 0, the wind-induced motion of a 
phase screen with the wind velocity V is to be taken into con¬ 
sideration since it causes extra averaging of phase fluctuations 
along a straight line of length VT. This effect is described by 
convolution of the function (3) with a rectangular function 


rect(z) = 


(VT) -1 , z< VT/2 
0, z > VT/2 


( 6 ) 


where z is measured along the wind direction which forms 
some angle s with the x-axis. 

In differential astrometry the position of the target object 
is measured with respect to some N reference stars. Let us 
consider a geometry of a stellar group in projection on the 
phase screen displayed in Fig. 1. The points B,, z = 1... N 
represent reference stars, the target star projection Bo is given 
by the vector r. Cartesian coordinates x,, y, of stars are equal 
to their standard coordinates scaled to the layer height h. In 
polar coordinates, the distribution of points B, is given by 
their linear separations from the target Sj = hp ,• where p ,■ is 
z-th star angular separation from the target, by positional an¬ 
gles i/jj measured with respect to the x-axis and by distances 
Sij — Js’j’ + s 2 j - 2 SjSj cos (if/j - i/jj) between the two points B, 
and B ; . 
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Fig. 1 . Geometric elements specifying a star group (circles) projection 
on a phase screen 


Differential displacement of star images AJ is obviously 
equal to the difference of £(r) function values taken at points 
Bo and B,-. In the trivial case of one reference star Bj, the dif¬ 
ferential displacement is A£(r) = J(r) - £(r + Si) = £(r) * Qfr) 
where 

G(r) = tf(r)-tf(r + Sl ) (7) 

is a difference operator. For multi-star reference frames the 
form of the function (7) is specified by distribution of reference 
stars which we include by using designation Q( r, s). 

Taking into account the above effects, we come to a com¬ 
mon expression 

A£(r) = * rect(r) * P(r)/cr * Q(r, s) (8) 

/7T or' 

inherent to the differential technique. 

Now it is useful to proceed to the analysis in the frequency 
domain where the power spectral density F((q, ip) of the quan¬ 
tity Aff is expressed as a product of four squared Fourier trans¬ 
forms F of functions convolved in (8) and the power spectrum 
(1). In the polar coordinates q , ip 

F((q, ip) = [d/(2;r)] 2 r 2 {Jr} r 2 {^} direct) 

xT 2 {Q}Ft(q) 

It is known that T 2 {d/dr'} — [2nqcos(ip - 6 )] 2 and 

!F 2 {rect} = sinc 2 [7ryr<j'cos(i/3 - s)] where 6 is the angle formed 
by the vector r' and x-axis and sinc(z) = sin (z)/z (e.g. Martin 
1987; Lazorenko 2002a). Angular parameters s (the wind 
direction) and 6 (direction of image motion measurements) 
present in the above expressions introduce anisotropic effects 
(Lazorenko 2002a), which makes the analysis rather compli¬ 
cated and requires a special discussion. However, there are 
some reasons permitting us to perform averaging of Eq.(9) over 
e and 6. With respect to 6 , this procedure is justified by the fact 
that, normally, the image motion variance is defined as its mean 
value measured along the x and y axes, hence F 2 {d/dr'} = 
2n 2 q 2 . Averaging over another angle s is related to variations 
of the wind direction in the vertical profile of turbulence and 
to its temporal variability during the exposure in a single layer. 


Numerical estimates (Lazorenko 2002b) show that even small 
±5° variations in the wind direction cause an effect very similar 
to a complete averaging over 0-27T range. The averaging leads 
to expression (Lazorenko 2002a) !F 2 {rect} = ( nVTqy x valid 
for long exposures T » D/V. Finally, note that due to the 
forthcoming integration in the frequency plane, the angular- 
dependent function T 2 \Q\ in Eq.(9) can be substituted by its 
average 

Q(q) = (2zr)-' f T 2 {Q(r,s))dip (10) 

Jo 

which is the filter induced by the reference star field. Then 
Eq.(9) takes the form 

A 2 c 

F((q ) = 2nVfq' 2 ~ PY ( qmq) (11) 

where Y(q) = < F 2 \P/cr } is the filter-function of the entrance 
pupil. Integration of the power density (11) over q yields the 
variance of differential image motion 

a^c r°° 

A 2 = --f Y(q)Q(q)q- 1 -P dq (12) 

v 1 Jo 


3. Image motion power spectrum attenuation 

A concept of image motion variance (12) reduction takes ad¬ 
vantage of the dependence of the power spectrum I\-(q) on fil¬ 
ter Y(q) and Q(q) shape which generally can be corrected. To 
explain the concept of the method consider a classic example of 
double star observations with a filled monopupil (5) when func¬ 
tions Y(q) and Q(q) are given by expressions (e.g. Lazorenko 
2002a) 

Y(q) = [2J | (n Dq) / (n Dq)] 2 
Q(q) = 2[l-Jo(27iq Sl )], 


where J m are Bessel functions of the order m. In this study we 
suggest some improved modifications of filters (13) for better 
filtration of turbulent phase distortions. For this purpose we re¬ 
quire that: 1) the low-pass filter Y(q) is fast decreasing at fre¬ 
quencies longer than D 1 , and 2) the high-pass filter Q(q) has 
very low response at short frequencies. A useful approximation 
of filter shapes is given by quasi-rectangular functions with flat 
peaks and drop-down segments: 


Y(q) = 


Q(q) = 


1, nDq<q Q 

E v (Dq/2)~ v , nDq > q Q 
H k (Sq) k , nSq<\ 

Qoc, nSq > 1 


(14) 


Here the terms v and k are the most important model param¬ 
eters since they determine the asymptotic behaviour of filters; 
E v , IIi and Q,x are quasi-constant terms, approximately inde¬ 
pendent of v and k; qo specifies the filter Y(q) nucleus width; 
S is the effective linear size of the stellar group projection onto 
the phase screen. 

In terms of this model, a very narrow field regime of ob¬ 
servations is defined by a nonoverlapping position of Q(q) and 
Y(q) nuclei: S D/qo. With an approximation go ~ 2 follow¬ 
ing from Eq.(18) it is simply S <K D/2. 
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Fig. 2. Combined filter Y(q)Q(q) response at D/S = 17.2 for the two 
configurations of Table 2: ”a” (upper curve) - double star, k = 2, v = 3; 
and ”g” (lower curve) - at k = 8, v = 9 


Classic filters (13) are described by parameters v = 3 and 

k = 2. 

A qualitative demonstration of the current concept of en¬ 
hanced image motion filtration requires using results from the 
next Sections 4-5. In Fig.2 we display the combined filtering 
effect of the product Y(q)Q(q) for a double star observation 
(k = 2, v = 3, upper curve) which should be compared with 
that expected for improved parameters v = 9 and k — 8 (lower 
curve). The functions plotted refer to a D - 100 m telescope, 
a single turbulent layer at h = 20 km and star reference groups 
”a” and ”g” of equal effective linear size S = 5.82 m (Table 2). 
The much lower filter response in the second case ensures the 
essential reduction of the differential image motion spectrum 
and A 2 . 

With the model (14), and the assumption of very narrow 
fields S «: D/2, the value of A 2 defined by the integral (12) is 
given by a two component sum 

A 2 = [c\(2S /D) k + ci(2S / D) V+P ](D/2) p n v ~ k+p (15) 

where ci and co are some slow functions of the model param¬ 
eters. Under condition k < v + p, the first component prevails 
resulting in a power law A 2 = (2 S / D) k D p while a dependence 
A 2 = (2 S/D) V+P D P is expected at k > v + p . This situa¬ 
tion is schematically reflected in Table 1 where the expected 
power laws for A 2 are shown. Diagonal elements of the table 
are marked with boxes and physically impossible k (odd) are 
omitted. 

Remembering that very narrow field 2S/D <sc 1 regime of 
observations is discussed, from Table 1 we find that the use of 
high v and k orders may result in a very strong reduction of A 2 
magnitude. Thus, while power laws A 2 ~ I) 2>p and A 1 ~ S 2 
are valid at k = 2 and v = 3, application of filters (14) with 
v - 9 and k = 8 leads to much stronger dependencies A 2 ~ 
D~* + p and A 2 ~ A 8 , owing to which about a 256-fold change 
in A 2 magnitude is expected for a 2-fold change in D or S. In 
comparison to double star observations, the expected decrease 
in A 2 is roughly (2 S /D) 6 . For the case shown in Fig.2, the gain 


Table 1. Asymptotic dependency of image motion variance A 2 on S 
and D 



! k 

V 

2 3 

4 

5 

6 

1 

8 

3 

S 2 D~ 2+P □ 

S 3+p D- 3 

- 

S 3+P D 3 

~ 

S 3+p D~ 3 

4 

S 2 D~ 2+p - 

S 4 D- 4+p 


S 4+p D- 4 

~ 

S 4+p D- 4 

5 

S 2 D~ 2+p - 

S 4 D- 4+p 

□ 

S 5+p D- 5 

~ 

S 5+p D~ 5 

6 

S 2 D~ 2+p - 

S 4 D- 4+p 


S 6 D -6+ p 


S 6+p D- 6 

7 

S 2 D 2+p - 

S 4 d- 4+p 

- 

S 6 D -6+ p 

□ 

S 1+p D- 1 

8 

S 2 D 2+p - 

S 4 D- 4+p 

- 

S 6 D -6+ p 


S s D-* +p 

9 

S 2 D 2+p - 

S 4 d- 4+p 

- 

S 6 D -6+ p 

- 

S g D-* +p 


is about (8.6) 6 « 4 • 10 5 . This estimate is rather approximate; 
more reliable assessments are given by Eq.-s(36) and (38). 

Table 1 shows that the optimal sets of k and v are concen¬ 
trated near the diagonal, so the best suppression of image mo¬ 
tion is expected at k ^ v. A further increase of one of these 
parameters, e.g. k at fixed v, leads to no improvement in power 
law though it may essentially affect the magnitude of A 2 . 


4. Apodization of the entrance pupil 

In this section we suggest an improvement of the filter Y(q) fre¬ 
quency response by applying a special apodizing mask in the 
pupil plane. It is known that a proper apodization of the en¬ 
trance pupil ensures faster decay of the Fourier transform T{P) 
of pupil function P(r) at high frequencies (Papoulis 1971). 
In application to astrometric measurements, apodization al¬ 
lows one to suppress high-frequency phase distortions which 
is equivalent to attaining a high v order of the function Y(q). 

It is known (Papoulis 1971) that when the function P(r) and 
its first n — 2 derivates become zero at the pupil’s edge r — D/2, 
then its Fourier transform decreases asymptotically as q~ n ~ 1 / 2 
at high q. The function Y(q) = !F 2 {P) in this case is approxi¬ 
mated by Eq.(14) with a parameter v = 2n + 1. To obtain high 
v orders, it is therefore necessary to use functions P(r) which 
are sufficiently flat at pupil edges. Also, apodization should en¬ 
sure sufficiently high light-transmission of the objective. For 
instance, definition of P(r) as a convolution of t he two func - 
tions (5), or P(r) = 7 G 1 [2 arccos(2 r/D) - (4r/D) y/l - 4r 2 /D 2 ] 
leads to Y(q) - [2J\{nDq/2) / {nDq/2)] A with v = 6 and an 
effective collective area (4) of the objective cr = nD z /l6. A 
light-transmission of the objective 

y = 4cr/(^D 2 ) (16) 

with this type of apodization is only y = 1/4. Better light- 
transmission is provided with the function 

P(r) = (1 - 4r 2 /D 2 y 1 , 


(17) 
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where n is an arbitrary positive integer. Functions Y(q) corre¬ 
spondent to (17) and their characteristic parameters (14) are 
given below: 


m = [2 cr=*-£, y = n\ 

E v = (nl) 2 /( 2tt 2 ' !+2 ), v = 2n + 1, q 0 = 


(18) 


An expression for effective filter width qo was found from the 
natural condition f f Y(q ) dq - n(qo/nD) 2 -, at moderate n qo ~ 

2 . 


The functions P(r) plotted for n — 1... 4 and corresponding 
to v — 3 ... 9 odd are shown in Fig.3. The filters Y(q) formed by 
this type of apodization are shown in Fig.4 and have asymptotes 

~ q- v - 



0 0.1 0.2 0.3 0.4 0.E 

r/D 

Fig.3. The pupil functions P(r) defined by Eq.(17) for n = 1.. .4 
(v = 3 ... 9) 



1 10 


q. P'] 

Fig. 4. Functions Y(q) correspondent to P(r) shown in Fig.3 

A negative, inevitable consequence of high v order imple¬ 
mentation is the decrease in light transmission y and so a weak 
light signal. Thus a certain trade-off exists between good filtra¬ 
tion of image motion and a high signal level that is necessary 
for image centroiding. 


5. Reference frames with a virtual symmetry 

Consider a vector quantity W with two components 

W x = A" 1 at(x o - Xj) = x 0 - A -1 Zti a <*i n 

W y = N - 1 ai(yo - yd = yo~ N~ l a,* ( ’ 

formed by a linear combination of measured (instant image 
motion included) differences of Cartesian coordinates of the 
target Jco, >’o and ;-th reference star Jcy,, i = 1,2..A scaled at 
the phase screen height h; a, are some weights that satisfy a 
normalizing condition 

N 

Yj a ‘ = N - ( 20 ) 

i= 1 

Introduction of the quantity (19) instead of normally non- 
weighted differential positions reduces the problem of image 
motion suppression to the qualification of conditions on a, 
which minimize the variance of W. A peculiarity of this study 
is the permission to use any, even negative, weights a,-. The vec¬ 
tor W that defines the differential position of the target object 
relative to the weighted center of reference group A" 1 £ 

A -1 2 a i}’i is the only quantity which can be measured with 
high precision. Therefore, the possibility to extract positional 
information from W is not obvious and will be discussed in 
Section 7. 

Besides atmospheric noise, the total error of observations 
depends also on the image-centroiding error component <t p i, in 
W caused by a Poisson photon noise in star images. Assuming 
that Dj is the variance of the ;-th reference star centroid position 
caused by photon noise, we find from Eq.(19) 

° 2 ph = N- 2 Y j a 2 D i (21) 

i= 1 

since a contribution from the usually bright target object is 
small. For Gaussian shaped images with r.m.s. width cr 0 = 
FWHM/2.36 and n, detected photons, an expression (Irwin 
1985) 

Di = a^/m (22) 

is valid for sufficiently bright images with a light signal exceed¬ 
ing photon noise. 

5.1. Functions Q(q ) of high k-orders 

Expressions for differential image motion effects derived in 
Sections 2 - 3 in application to the double star observations 
can be easily extended to the case of a measured value W. 
Image motion for the quantity (19) is now expressed through 
the weighted differences of the function C values in points Bo 
and B,-: Ac = A -1 £, a;[(T(Bo) ~ GAj]- The difference opera¬ 
tor (7) corresponding to and the geometry shown in Fig.l is 
then 

Q(r, s) - N- 1 ai [6( r) - <S(r + s,)] (23) 
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A Fourier transform of this expression with a subsequent inte¬ 
gral averaging ( 10 ) yields a function 


N 

Q(q ) = N" 2 Y «;«;[! - 2J 0 (2nqSi) + J 0 (2nqSij)] (24) 

U= 1 


corresponding to a new definition of the measured quantity. 
Analysis of Eq.(24) permits us to formulate the conditions nec¬ 
essary to increase the Q(q) function order. For this purpose we 
expand the Bessel functions in (24) into power series of q , de¬ 
riving the approximation at low frequencies: 


Q(q) = wZ 


m= 1 (m!) : 


;j-0 *q) E ?/=1 a/flf 


(25) 


x(s : 


2m . „2m 


' #> 


The order A' of the Oft/) function is determined by the first non¬ 
zero term in the expansion; normally, k = 2 if no special efforts 
are made. 

With Sij fixed, one may always choose a set of co¬ 
efficients a, turning to zero the first or few first sum val¬ 
ues Ef,= i a i a j( s f n + s )'“ ~ s fj ”) which are coefficients at q 2 ’". 
With a sufficiently large number of reference stars, all the first 
q 2 , q 4 ... q k 2 expansion terms in (25) could be eliminated by 
applying conditions 


2 E atajis 2 + s 2 - s 2 j) = 0, 
E E aicijis 4 + s 4 - s 4 j) = 0, 


E E mA^T 2 + s k r 2 - sp 2 ) = o 

where k is some even integer. Under conditions (26), the power 
of the first non-zero term in the expansion (25) can be increased 
from normal k = 2 to some higher k > 4. Eq.-s(26) is a system 
with unknowns ar, it includes only a single first line at k - 4, 
the two first lines at k — 6 , and so on. The system (26) can be 
converted to a simpler form. Thus, passing to Cartesian coor¬ 
dinates Xj - vo = Si cos ijfj, yi - yo — Sj sin if/i and taking into 
consideration Eq.(20) one finds that each sum of sums in Eq.- 
s(26) takes the form of products of one-dimensional sums. For 
instance, the first equation leading to k = 4 transformes to 

Y X a ‘ a J^ + s2 J~ + M y) = 0 (27) 

where M x = E a,(x,—xo) and M y = E afyi—yo) are the first co¬ 
ordinate moments. The next equation, for k = 6, E E a i a j(s 4 + 
s 4 j - Sp = 2[E cij(xi - xo) 2 + E«iCy i - yo ) 2 ] 2 + 4[E«;(*i - 
•*o) 2 ] 2 + 4[E afyi - y 0 ) 2 ] 2 + 8[E a^x, - x 0 )(y,- - y 0 )] 2 = 0 con¬ 
tains quadratic cross-moments. Direct computations show that 
the order of cross-moments is incremented by 1 when passing 
to each next k order. The quadratic form of the equations sug¬ 
gests that all weighted cross moments of reference star coordi¬ 
nates are zero. An equivalent form of system (26) is therefore 


E a i = N, 

E a,(x, - X 0 ) = E a i(yi - Fo) = o, 

E flj(Xj - x 0 ) 2 = E a-iiXi - x 0 )(y; - y 0 ) = 

= E a i(}’i - yo ) 2 = o> 

E a,(x,- - x 0 )i _1 = E «;(•*'; - -*'o)^“ 2 (y/ - yo) = 

E a,(x,' - X 0 )(y; - yo)^ -2 = E afyi - yo)^~' = 0 


(28) 


which reveals the modal structure of the filter Q(q). At k = 2, 
when none of conditions (26) are fulfilled, the system is limited 
by a single first line normalizing equation; at k = 4 it includes 
the first two lines (three equations), and so on; the total number 
of equations is k(k + 2)/8. In a compact form, the system (28) 
with unknowns a, is written as 


E a i = N, 

E aj(xi - x 0 ) a (yj - yof = 0, a +/3 = 1... 5 - 1, 


(29) 


where a and /3 are positive integers. The use of weights a,- sat¬ 
isfying Eq.-s(29) results in elimination of A/2 - 1 first terms in 
the expansion (25) which becomes 


0( a \ - _L y°° a im V iv 

Siyq> - N 2 2-‘m=kl 2 ( m |)2 H 

x(s 2 


i,j= 1 a i a j 


(30) 


Plm „ 2 m 


„ 2 m^ 

O' ; 


It should be noted that in the case of k > 6 when quadratic 
moments of x and y are involved in Eq.-s(29), the solution must 
incorporate negative a, values, which is rather unusual for a 
common technique of astrometric reductions. 

Application of various sets of a,- changes the geometric 
properties of a reference group. Thus, implementation of a, 
of the k = 4 order for which condition E «/(x,- - xo) = 
E afyi - yo) = 0 is fulfilled transforms an arbitrary reference 
group into its virtual equivalent with an ideal symmetric struc¬ 
ture centered at the point xo, yo- With respect to the image mo¬ 
tion statistics, and the filter Q(q) asymptotic behaviour in par¬ 
ticular, both groups become indistinguishable providing they 
are of equal effective size. A further increase of k results in 
enhanced improvement of reference group filtering properties 
which can reveal a power dependency stronger than Q(q ) ~ q 4 , 
a case impossible with a simple geometric symmetry. 

Implementation of high k orders is limited by the number 
N of reference stars available. The minimum N value required 
to achieve some k order and to find solutions a of the system 
(29) is 


_ j k(k + 2)/8, for 2-D distribution of stars 
Nmin ~{k/2, for 1 -D, in-line distribution (31 ’ 

Here we admitted that N mm , generally, depends on the type of 
star distribution in the field. Though exactly linear configura¬ 
tions do not exist, they serve well for illustrative purposes. 

At N > N m i„, a redundancy of the system (29) allows us to 
set a useful (even in the case of k — 2) condition 

y a 2 Dj = min (32) 

which reduces the centroiding error (21). A consistent solution 
of (29) and (32) is found with a standard Lagrangian method 
of undetermined coefficients. Such a solution is optimal with 
respect to both atmospheric and photon centroiding errors. 


5.2. Tutorial linear configurations 

Some tutorial examples of linear (along the x-axis) configu¬ 
rations are given in Table 2 which contains the conditional 
name of configurations, x-coordinates of stars given with ref¬ 
erence to the target and expressed in seconds of arc, a, values. 
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a 2 = N 1 Y aj and the principal q k term of the Q(q) function 
expansion, x,- are scaled to equalize the effective angular size of 
each group to p = 1'. For a turbulent layer at h = 20 km, this 
angle corresponds to an effective linear size S — ph — 5.82 m. 
The last quantity is convenient to define as 


S = 


N~ 


[ ^ didjis) 




i/k 


(33) 


which relates S to the magnitude of the first non-zero term in 
the expansion (25). A convenience of this definition is that se¬ 
ries expansion (30) (its leading term) is now simplified to 

Q(q) = (nqS'f/[(k/ 2)!] 2 (34) 


Note that effective sizes S and p depend very weakly on the 
peculiarities of the star distribution in the frame and approxi¬ 
mately are equal or slightly exceed its radius (the largest sepa¬ 
ration target-reference star). 



Fig. 5. Functions Q(q) for stellar groups ”a, c, d, g” (Table 2) with 
k = 2, 4, 6 and 8 symmetry order; D/S = 17.2 


Table 2. Linear configurations of effective angular size p = 1' and 
k-order virtual symmetry 


Name N 

k 

X; 

a./ 

a 2 

Q(q) 

a 

1 

2 

42" 

i 

i 

f 

b 

2 

4 

-38, 38 

i, i 

i 

q 4 

c 

2 

4 

27, 54 

4, -2 

10 

q 4 

d 

3 

6 

-29, 29, 58 

1,3,-1 

3.67 q 6 

e 

3 

4 

-32, 32, 64 

12/7,6/7,3/7 

9 

q 4 

f 

4 

4 

-31 V2, -31, 







31,31 V2 

1, 1, 1, 1 

1 

q 4 

g 

4 

8 

-30 a/2, -30, 







30, 30 VI 

-2, 4, 4, -2 

10 

q* 

h 

6 

12 

-56,-37.4,-18.7, 

0.3,-1.8,4.5, 






18.7,37.4,56 

4.5,-1.8,0.3 

7.9 

q 12 


Table 2 shows that even strongly asymmetric star groups 
”c” and ”d” are subject to high order symmetrization. This is 
achieved, however, at the expense of applying large a, owing to 
which the a 2 value is greater than that for symmetric groups 
(compare cases ”b” and ”c”). Thus, application of high or¬ 
der symmetry for asymmetric star groups reduces atmospheric 
noise but causes a rise of the centroiding error (21). The same 
effect is observed even for symmetric distributions ”g” and ”h” 
when very high k is used. 

The Q(q ) plots shown in Fig.5 for a few configurations of 
Table 2 emphasize a difference in the function form at low 
frequencies which for both symmetric and non-symmetric star 
distributions is determined only by the index k. Asymmetry of 
star groups ”c” and ”d” leads to the increase of amplitude Q m 
°f Q(q) at high frequencies q > 1 /S, which, due to relation 
a 2 /N = Qoo - 1 following from Eq(24), is indicative of the cr 2 ph 
increase. 

An example of 2-D stellar group symmetrization is given 
in Table 3 for a sample of N — 14 stars of the open cluster 
NGC 2420 observed with the 5-m Palomar telescope (Pravdo 
& Shaklan 1996). The table contains star coordinates x,, y, with 


reference to the target star, weights a,- and d 2 - Y d 2 /N for 
various orders up to k = 8. Weights a,- were computed using 
conditions (32) with D, = const. The extremally small angle 
p = 3.7" found at k — 2 is by no means due to a special selec¬ 
tion of reference stars (no selection was applied) and reflects 
the effect of field averaging making the first moments M x and 
M y small for large N. In the limit of /V —> °o, discussed in 
Section 6, any k — 2 order group is moved up to the 4-order. 
Note also that the value of a 2 is rapidly increasing with k , so 
does the variance (21). 


Table 3. Weights a, symmetrizing a stellar group (Pravdo & Shaklan 
1996) to k = 2... 8 orders 




1 k 

X; 

yi 

2 

4 

6 

8 

46.0" 

19.4" 

1.0 

0.914 

-1.355 

-0.023 

33.5 

17.1 

1.0 

0.908 

0.340 

-1.120 

19.1 

18.7 

1.0 

0.845 

0.739 

0.680 

9.5 

16.1 

1.0 

0.851 

1.792 

0.308 

-3.2 

5.7 

1.0 

0.951 

3.414 

3.300 

17.1 

3.4 

1.0 

1.039 

3.201 

4.300 

OO 

OO 

-11.4 

1.0 

1.209 

1.787 

7.102 

14.6 

-15.7 

1.0 

1.282 

0.624 

-0.265 

19.2 

-16.1 

1.0 

1.300 

0.465 

-2.793 

-24.1 

1.1 

1.0 

0.950 

2.981 

2.157 

-32.5 

-4.7 

1.0 

1.002 

1.940 

2.362 

-41.3 

-8.6 

1.0 

1.028 

0.572 

-1.593 

-36.3 

-16.4 

1.0 

1.144 

-1.407 

-0.138 

-45.5 

24.9 

1.0 

0.578 

-1.092 

-0.278 

a 2 

1.0 

1.03 

3.40 

7.32 

p 

3.7" 

45.7" 

35.2" 

42.2" 
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5.3. Expressions for the variance of image motion 

Analytic expressions for A 2 in a limiting case of very long ex¬ 
posures and narrow fields T » D/V » S/V can be derived 
from Eq.(12) where Y(q) and Q(q) are given by Eq.-s(18) and 
(24): 


A 2 


= a^W!) 2 

n 2 vt 


OO 9 

[•WW 
J (nDq) v ~ l 

0 


E E«;0/[ 1- 


-2J 0 (27rqSi) + J 0 (2nqsij)]q 1 p dq 


(35) 


An approach to integration depends on the ratio between the k 
and v values. When k > v+p, the integral converges at q - 0 so 
approximation (14) for Y(q) is valid. Integrating by parts and 
taking advantage of the condition (26) results in 

,ty(^i!)>v f T(*±i)T(^), V- odd 

An^vTD^rqe^) \ T(f )T(^), v- even (36) 

k> v + p 


Here S — |/V“ 2 X E a,a ; [2,v, p+v - .v ! ■/ ,+v ]| l ^ ^,+ ' ,| is a modified 
field size defined similarly to (33) but with fractional pow¬ 
ers. Since s, = sp, /p, from (36) we find a power dependency 
A 2 ~ S v+P /D v which is in agreement with Eq.(15) and Table 
1 testifying that an increase of k over v + p does not affect the 
power laws. Another interesting point is the dependency of A 
on k which, however, is not clear from (26) since it contains 
the S term related to k via weights a,-. A rough expression for 
A with no S term can be obtained based on the next Sections 
results for non-discrete star fields. Using approximation (49) 
for Q(q) and recomputing the integral (12) with Y(q) given by 
(14), yields 


2 _ iV y - 1 (ty4) t - v -^!)V + p 

VT[(k/2)\] 2 (k-v-p)D* 


,k > v + 1 


(37) 


which is a decreasing function of k at v fixed. The use of high 
k orders, thus, always reduces atmospheric noise. 

At k < v + p, the integral (35) is calculated with approxi¬ 
mation (34) for Q(q). Direct integration yields 


2 _ 

vnat/2)!] 2 ^D‘-rr(v+£^)r(ti£ik±) (38) 

k < v + p 


and is also consistent with Table 1. It should be stressed that, 
unlike the previous case of k > v + p, an increase of v over 
k gives rise to A 2 . The deterioration of the results is related to 
the expansion of the filter Y(q) nucleus width, clearly seen in 
Fig-4. 


6.1. Approximate expressions for at and cr ph 

We assume equal brightness of stars. Assuming also that refer¬ 
ence stars are evenly distributed around the target in a circle of 
radius R with a spatial density N/(nR 2 ), we introduce, instead 
of weights a,, a weighting function ci(r) with a radial symmetry 
which satisfies the normalizing condition f () a(r)dxdy = 
N equivalent to (20). For a /.-order function Q{q), an integral 
analogue of the system of equations (29) with unknown func¬ 
tion a{r) is 

f*a(r)rdr = R 2 / 2, k = 4 

f 0 R a(r)r a dr = 0, a = 1,3 ... ^,k = 8,12... 

The above system is valid for k multiples of 4 only, since for 
any symmetric distributions Eq.-s (29) with coordinate cross¬ 
moments of odd powers are satisfied automatically; for this rea¬ 
son orders k - 2,6,10... do not exist, they are moved up into 
the next higher order. Considering only a polynomial class of 
solutions for a(r) of the form a(r ) = bo + b 2 r 2 + ... + bprP, 
where / 3 = k/2 - 2, and performing integration (39), we come 
to a linear system of/3/2+1 equations with respect to unknowns 
b 2 f 


Z 


j =o 


bijR 21 
2(i + 7+1) 


R 2 , i = 0 

0, i = 1,2.. .yS/2 


The solution found by computer simulation is 


„( r \ - k yk/ 4/ iy+1 (ft/4-i+l)! 

4 (Jfc/4-i)!i'!(i-l)! 

k = 4,8...; r<R 



(40) 


(41) 



r/R 

Fig. 6. Weighting functions (42) for circular star distributions at k = 4, 
8, 12; the asymptotic form (43) of a(r ) at k = 12 


6. Virtual symmetry for dense reference frames 

Approximate characteristics of reference groups with several 
stars are easily found in a limit of infinite N when the star distri¬ 
bution becomes continuous. All estimates are found especially 
easily since discrete summations are substituted by integrals, 
and individual features of star distribution in the field become 
unimportant. 


For a few first k. the weighting functions 

(1, k = 4 

a(r) = \ 2(2 - 3 r 2 /R 2 ), k = 8 (42) 

( 3(3 - 12 r 2 /R 2 + 10 r 4 /R 4 ), k= 12 

are shown in Fig.6. Except for the case of k = 4 dis¬ 
cussed by Lindegren (1980) and when air) =const, the plots 
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of a(r ) incorporate some oscillating circular zones coming 
from the field outer borders. The number of these zones 
and the function’s peak increases with k while the ampli¬ 
tude of oscillations decreases. These features are easily ex¬ 
plained since Eq.(41) at k —» oo is reduced to air) = 
(. k/4 ) 2 Zfio~'(-l)'[kr/(4/?)] 2 7[;!(i + 1)!] which is a truncated 
power expansion of the Airy-type oscillating function 


« 0 ) = 


(k\ 2 2Mg) 
14 ) kr/(2R) 


(43) 


For comparison, the asymptotic and exact forms of a(r ) for k = 
12 are shown in Fig.6. 

Above we assumed equal brightness of field stars. Fet us 
find expressions for a weighting function a(r, D,) that approx¬ 
imates weights di for descrete very dense distributions while 
accounting for magnitude-dependent centroiding variances Z), 
(22). In this case the condition (32) leads to the inverse propor¬ 
tion «, ~ 1/Di, which follows also from a least-squares princi¬ 
ple. Assuming that spatial and brightness distribution of stars 
in the field are uncorrelated and using a(r) as approximation 
for a, at Z), =const, we obtain 

a(r,Dj) a(r)^ (44) 

where D\ = N / E/Ii ZU 1 is the mean (effective) variance of one 
reference star. The total centroiding variance of all reference 
field than is Z9 e ff = D\/N = 1 /2Z>7* . Taking Eq.(22) into 
consideration, we find 

Oeff = cr ( 2 / ^ m. (45) 

Incorporating fainter stars, a value of Z) e ff always decreases 
(improves) but rapidly approaches some limit. 


20°, l = 0°. The field with stars to V=23 mag was symmetrized 
to the k = 12 order; a distribution of computed a, values versus 
distance from the field center is shown in Fig.7. In this exam¬ 
ple an accumulated value of Z9 e ff is equivalent to that produced 
by 120.3 stars of 15 mag; an effective mean centroiding vari¬ 
ance of one reference star Z>i corresponds to the star of V= 18.5 
mag. A large scatter of computed a, values seen in Fig.7 is nat¬ 
ural since a ; - 1 /£),-; this relation is exact at k = 2 and becomes 
rough at large k. Usually, |a,-| » 1 for bright and |a,-| <sc 1 for 
faint images. Sampled values of a,D,/D\, on the contrary, show 
small scatter and a strong concentration around the a(r) func¬ 
tion plot (Fig.8). Approximate values of a, thus can be found 
from Eq.(44) for any k star. 



star radial distance, arcmin 

Fig. 8. ciiDj/D i values for random data shown in Fig.7; solid line - 
approximation by the function a(r) plotted for k = 12 



star radial distance, arcmin 


An increase of the air) amplitude in the field center for high 
k leads to an increase of the total centroiding error <r 2 ph . This 
effect is easily estimated since, with above the assumptions, 
the magnitude-related D, and coordinate-related a 1 terms in 
Eq.(21) are statistically independent. Therefore Eq.(21) trans¬ 
forms to cr 2 ph = D^a 2 where a 2 — (zrZ? 2 ) -1 f f a 2 (r)dxdy is 
an averaged value of a(r). Direct integration of Eq.(42) yields 
a 2 = (k/4) 2 whence 

cr 2 ph = D eS (k/ 4) 2 (46) 

The last equation can be given in terms of star image param¬ 
eters. Taking into consideration Eq.(22) and allowing for light 
transmission y of the apodized objective Eq.-s(18), we obtain 


Fig.7. An example of a, coefficients symmetrizing a random 5' field 
with 3170 stars to the k = 12 order; dot sizes are proportional to the 
star brightness (12 to 23 mag) 


tr ph — 


FWHM 
2.36 VJT); 


(k/4) 



(47) 


The above considerations are illustrated by a numerical 
simulation of a random 5' field (radius) with 3170 stars. The 
magnitude and sky star distribution corresponds to the Galaxy 
model (Bahcall & Soneira 1980) and galactic coordinates b = 


The photon noise thus depends on FWHM, the total light of 
reference stars E n h k and v. The use of very faint stars as a 
reference clearly does not lead to improvement in cr ph because 
their contribution to E f U is negligible; the use of very high k 
orders degrades cry,. 
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6.2. Expressions for Q(q) and S 

The asymptotic form of the function Q(q) for stars of equal 
brightness is found from a limiting (N —* °° ) form of the func¬ 
tion (23) which is Q( r) = <5(r) - P(r)ci(r)/(nR 2 ). From this it 
follows that Q(q) = T 2 {Q) = [1 - 2 R 2 f* a(r)J 0 (2nrq)rdr] 2 , 
and a direct integration with a(r ) given by Eq.(42) yields 

' [1 - 2J\ (2nRq) / (2nRq)] 2 , k = 4 

(i 4J l (2nRq) _ TAJiOfRqp 2 i. _ q 

L 1 + (2 nRq) (2 nRq? \ ’ K ~ 8 


that gives the expectation of effective field size for a random 
star sample taken in a circle of a radius R\ both S exp and R 
are given in arcminutes. The model (51) is consistent with (50) 
and extends this dependency to odd kj2. The quantities^ and 
ft valid for both high and low star density are given in Table 4. 
At polar regions, the values of these parameters are higher, thus 
slightly larger S and A are expected at equal R. 

Table 4. Coefficients of Eq.(51) for different star density 


Q(q) = 


, _ 6J l (2nRq) 96J 2 QxRq) 
(2nRq) (2 nRq) 2 

480J 3 (27rfej) t 2 
(2nRq? \ ’ 


0, q < k/(4nR) 
1, q > k/{47iR) 



| near-equatorial zone 

| Galactic pole 

k 

X 

R 

X 

R 

2 

0.275 

0.414 

0.848 

0.537 

6 

0.712 

0.817 

1.000 

0.830 

10 

0.857 

0.893 

0.997 

0.922 


X = fJ = 1 for k = 4, 8,12 (k/2 even) 


An expression for k —» oo derived with a Hankel transform 
of Eq.(43) shows that an asymptotic form of Q(q) is an opaque 
circle with radius proportional to k. Expansion of Eq.-s(48) into 
power series of q yields an expression 


Q(q) = 


q < k/(4nR ) 
1, q > k/(4nR) 


(49) 


valid for k - 4,8 _Comparing this expression with a defini¬ 

tion (34) for effective size S , we come to a very simple relation 


For k/2 odd the value of S is a function of odd-order coor¬ 
dinate cross-moments averaged over the field. Result of averag¬ 
ing depends largely on the contribution from a few bright stars 
with large a, and so is not returned to zero as it is expected for 
star fields with constant star brightness. Not performing a de¬ 
tailed study, we assumed that in the spectral domain this effect 
is described by the empiric model 


S = R, k/2 even 


(50) 



Fig. 9. Effective field size S r „d versus radius R of random samples at 
high star density; k is the reference group symmetry order; lines are 
data fits (51) 

To test the validity of Eq.(50), we performed numeri¬ 
cal simulations that assumed the Galaxy model by Bahcall 
& Soneira (1980). Random star samples taken in the near- 
equatorial zone (b = 20°, l - 0°) to V=23 mag were sym¬ 
metrized, with condition (32), to k — 2,4.. 12 order and esti¬ 
mates S rn d of 5 for each random sample were computed using 
Eq.(33). Random sizes S rn d presented in Fig.9 are fitted with 
the function 

(51) 


Qiq) = Q k+ 2 (q ) + {SIR) k [J 2 kl2 (2nRq)+ 

+2 2 J 2 m (2nRq)], k/2 odd (52) 

m=k/2+l 

with correct asymptotic properties. Here Qk+iiq ) is a k + 2 or¬ 
der function (48); S is given either by Eq.(51) to produce the 
mathematical expectation of the filter or is a sampled field size 
leading to a ’’sampled” filter. The second component in the ex¬ 
pression is caused by incomplete averaging of odd-order co¬ 
ordinate moments and decreases with R since (S/R) < 1. For 
fields with stars of equal brightness, this component vanishes 
yielding Q{q) = Qk+ 2 (q) and thus increasing the field symme¬ 
try order by +2. It is easy to find that Eq.(52) at short q follows 
approximation (34), critically important for computation of A. 
Direct calculations of A for random star fields in Sect.9 had 
proved the validity of the model (52). 

Finally let us derive laws for A as a function of refer¬ 
ence star number N in the field, and its radius. As follows 
from Eq.(31), the field symmetry order not can exceed k max = 
V8 N + 1 - 1. Then, with some optional k < k m . dx and optimal 
apodization parameter v-kovv-k+l (Table 1), we find that 
for a turbulent layer at a height h 




hR 


3400 D 


R 


Mm-D d 2/3 


(53) 


where R is in arcminutes, D and h in meters. Consider that at 
high star density, a R = V field is expected to contain about 
100 stars to V=23 mag; using all of them as reference makes k 
orders to k max ~ 25-30 quite feasible. A gain in A for narrow 
hR < 3400/) fields is therefore huge. 


S exp = x Rtl 
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7. Plate reduction 

An approach to plate reduction based on the current technique 
should take into account the fact that effective filtration of at¬ 
mospheric noise occurs only in a two-component vector W. As 
opposed to the common technique based on plate constant de¬ 
termination that provides a global fitting of measured to model 
data at any point on the plate, the present technique gives a local 
solution valid directly at the point xo, yo. When multiple targets 
are considered, each one should be processed with its own set 
of weights cij centered on the particular target. Below we ex¬ 
plain in detail the feasibility of recovering positional informa¬ 
tion from the quantity IT, taking into account the two possible 
goals of observations: 

- determination of the target object position with reference 
to field stars, using their precise positions (a classic problem of 
astrometry), and 

- determination of the target proper motion with reference 
to field stars (positions unknown) from observations made in 
the two epochs. 

7. 1. Compensation of low-order reduction model 
terms in W 

Consider polynomial expansion 

Xi = Xj + A q + A\'°(xj - x 0 ) + A^tyi - y 0 ) 

+A 2 2 °(xi - x 0 ) 2 + A l 2 ’ l (xj - x 0 )Of -y 0 ) 

+A®'"(y; - yo ) 2 + ... = Xj + Ao 
Jfc/ 2 -l m „ (54) 

+ E E - x 0 ) a (yi - yof 

m— 1 a,(3—0, 
a+(3—m 

+<Rki 2 {xi - xo ,y; - yo) 

of star measured coordinates xj, y, over coordinate cross¬ 
moments m = 1,2 ... A/2 - 1 of actual coordinates x,-, y,- not 
distorted by the atmosphere; the equation for y has a similar 
structure and all coordinates are considered to be scaled to the 
turbulent layer height h. Expansion (54) is essentially a reduc¬ 
tion model used for fitting measured to standard coordinates. 
The first few modes of the expansion with amplitudes A"f rep¬ 
resent common geometric terms: zero-point (mode m = 0 ) and 
scale (m = 1 ), also, they can be associated with the influence of 
the classic optical aberrations: tilt (m - 0 ), defocus and astig¬ 
matism (m = 1 ), etc. 

Each expansion term A ‘if may include also a stochastic 
component of image motion which, due to a similar structure 
of Eq.-s (54) and (29) is associated with a corresponding modal 
term of the Q(q ) function. The model is truncated by some 
k/ 2-1 order; the sum of high m > k/2 modes is denoted as 
2- 

An expression for the W x quantity is found by subtracting 
from (54) an identity xo = Ao + xo valid at the point x, = xo, 
y,- = yo and performing a summation with a,: 

W x = ^ E ai(x o - 

- -j/ E ai<Rki 2 {xi - xo,y, - yo) 

Here all A], A 2 .. .A/./ 2-1 low-order components vanish due to 
conditions (29). The use of weights a, based on the true dif¬ 
ferential position of stars xo - x,-, yo - y; therefore allows us to 


(55) 


form a linear combination (55) of measured coordinates which 
is insensitive to low-order terms of expansion (54), irrespective 
of whether they are constant or stochastic ones. It is very im¬ 
portant that W x does not depend on any changes of aberrations 
and temporal variations of image motion of low orders. The 
mathematical expectation of W x thus is zero and its variance 
depends on high k/2, k/2 + 1 ... modes uncompensated in the 
summation. 

7.2. Determination of target position 

Formulation of this particular problem implies that the precise 
position x 0 bj, y 0 bj of the target object is unknown and xo, yo is a 
preliminary object’s position used for computation of weights 
cij. The true target position x 0 bj = xo + Sx (6x is a correction 
to xo) is related, as it follows from (54), to a measured position 

Xobj - xq+Sx+Aq+A j ,0 <3x4 -A ‘I' 1 Sy+A 2,0 Sx 2 + Then Eq.(55) 

compiled for x 0 bj instead of xo becomes 


W x — A 1 E; Qiixobj - Xj) = 6x + A['°<5x + A"’ L 6y 

+ . . . - A -1 El 

Eq.(56) with a similar expression for y and conditions (29) for 
a, form a system with unknowns 6x, Sy. In the first approxima¬ 
tion, 6x = A -1 E; aifxobj - xi) yielding 

Xobj ~ Xo + Xobj — A ^ a,x, + A ^ ^ llSf^k /2 (57) 

For large Sx\ the system is solved by iterations, assuming A\ — 
0 ; weights a, are recomputed after each refinement of Sx, Sy 
and with a following shift of the point Xo, yo to a new position. 
Iterations converge very fast as A 1 «: 1. 

The above procedure allows us to obtain precise positions, 
for example, of extragalactic radio sources in the system of 
some high-accurate (future space mission) reference catalogue. 
Position of reference stars are used both for determination of x, 
y frame origin (reference group zero point) and for computa¬ 
tion of cij. Errors in reference star positions affect computed a, 
values thus implicitly causing a r.m.s. noise A a in x 0 bj positions 
(57). It is very important to know how large the A a component 
can be. 

To study this effect, assume that Sx,, Sy ; are random and 
uncorrelated coordinate errors of the ;-th star. Then the use of 
Xj + Sxj, }j + Syi data yields biased estimates at + Sa which 
results in inaccurate compensation of atmospheric error due 
to violation of conditions (28). In particular, the first moment 
M x = E(«; + Sai)(Xj - xo) is now not zero. Its value is easily 
found since the second equation in (28) used for computations 
of weights now takes the form J](ai+Sai)(xi-xo+Sxi) = 0 . Here 
the second-order Sci/Sx, terms can be discarded yielding M x = 
- Yj( a ii+Saj)Sxj ~ - E djSxj. Therefore, the q 2 term of the Q(q) 
function equal to 2(M 2 + M 2 ) = 2[(E ciiSxi) 2 + (E rijSyj) 2 ] in 
the left part of Eq.(27) becomes not zero. With respect to im¬ 
age motion statistics, its influence is equal to that caused by a 
q 2 term of the Q(q) function for a stellar group with effective 
coordinates x' = a t Sxj, yf = ci/Syi and unit weights. This group, 
in turn, can be substituted by a double star having a Q(q) func¬ 
tion (13) with distance parameter = V2 E fl ; cr / where <x, 
is the mean coordinate error of z-th star defined by equation 
Sx 2 + Sy 2 = 2a 2 . 


lO.l, 


(56) 
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Thus, inaccuracy in reference star coordinates results in 
the occurrence of an extra stochastic component whose be¬ 
haviour is similar to the differential image motion of an ’’equiv¬ 
alent” double star with a very small separation s\. Assuming 
equal star brightness, and taking into account that a 2 ~ ( k/A ) 2 
(Sect. 6 ), we derive s 2 « k 2 N<r 2 /8 where & 2 is the mean value 
of 07 . 

Consider, for instance, observations performed with a D = 
100 m telescope, T = 10 min, k = 8 and v = 9 in R = V field 
located close to the galactic plane (b = 20 °, l = 0 °), and with 
GAIA space catalogue positions as reference. In this sky area, 
the expected number of stars to the GAIA’s limit V=20 mag is 
N « 35. Assuming & « 20 /.(as as the average (Perryman et 
all. 2001) leads to si = 0.34 mas. Now, scaling Table 6 data 
with a dependency A 2 ~ S 2 valid for k = 2 (38), we find that 
for a double star system with S = si V2 = 0.48 mas, the noise 
A a is only 6 • 10 4 //as. It is easy to ascertain that for any other 
parameters of observations A a is always much smaller than the 
zero point error &/ V/V. 


7.3. Determination of proper motions 


The most important applications for astrometry with large tele¬ 
scopes, of course, are related to proper motion (exoplanet 
search) works. In this case the plate reduction technique is 
aimed at measurements of small displacements Ap of scien¬ 
tific objects which occur in the time interval between the two 
epochs T\ and 74. A frame of reference is given by the mea¬ 
sured positions of field stars x,, y,-. 

Extracting of high-accurate proper motion data is ensured 
by the next procedure. 

Observations of epoch T\ are used for computation of 
weights di(Ti) which correspond to some k symmetry order. 
A virtual center of the frame is fixed at the observed target 
position xo(7’i), yo(74) of the object location (sometimes, to 
avoid confusion, the epoch of a quantity measured is given in 
parentheses). Since true positions x,-, y, of stars are unavailable, 
computations are performed with the modified system 


Hi a; = N, 

Hi a-i{Xi - x 0 )“(y, - yof = 0 , a + p = 1 ... § - 1 


(58) 


using star positions shifted by the image motion. Solutions 
11,(7 1 ) therefore differ from unbiased solutions a,( 7 )) of system 
(29). Then, with a,(7’i) and measured star coordinates xo(74), 
xj(74) the epoch 74, we form 

W, = N- 1 Yj di(7’i)[x 0 (74) - x,(74)] (59) 


To estimate M x and M y magnitudes, assume temporarily 
that 74 refers to another time moment of the first epoch so that 
A/j - 0, and consider a model inversed to (54) and of a similar 
structure: x,- = x t + A 0 + A]’°(x; - x 0 ) + A®’ 1 ^,- -y 0 ) + • ■ • + % 2 - 
Assuming the model to be written for the moment 74, sub¬ 
tracting an expression xo = xy + ,4 (l valid for the target object, 
performing summation with weights a,{T 1 ) and taking into ac¬ 
count conditions (58), we find M x = H a,(7’i)[x I (7’i)-Xo(7’i)] = 
Hfl/(7’i)%/ 2 (7’i). A similar expression, of course, is valid for 
the moment 74 with ’R/q 2 related to 74. The quantity M x thus is 
a stochastic variable whose instantaneous value depends on a 
particular set of a,-. The mathematical expectation of M x as zero 
is reached at a, = a,-, its variance depends on the variance of the 
% 2 term and is therefore equal to the x-component of A 2 ex¬ 
pected at some current k and v. The average values of M 2 + M 2 
are thus of the order of A 2 . 

It follows that uncompensated extra image motion caused 
by a small q 2 component of Q(q) can be approximated (see 
a similar discussion in Section 7.2) by the image motion in 
the ’’equivalent” double star system with a separation ,v 1 = 

J M 2 + M 2 = A 2 . With A given in Table 6 one can evaluate 
si for specific parameters of observations; thus for a 10 m tele¬ 
scope we find si <0.1 mas if k > 4. Image motion induced 
in this double star system is very weak and normally can be 
disregarded. 

The use of measured star positions for computation of 
weights thus maintains the high accuracy of the method. 

In the above analysis 74 was related to the first epoch to null 
a proper motion effect. Putting 74 in the second epoch presents 
some problems since now xo(74) = xo(7’i) + A \x x . Using a re¬ 
versed model of measured to standard coordinates transform, 
we obtain the relation xo(74) = xo(74) + A\jl x + Aq + A['°A// V + 
A'j 1 ’ 1 A// v +... where A fi x and A/r v are the measured x and y com¬ 
ponents of A/j. Assuming for field stars x,(74) = x,(T 1 ), taking 
into account above expression and Eq.(58), we find 

AT‘_H ai(Ti)[xo(T 2 )- x,(7’ 2 )] 

= A// v + A]' 0 A// V + Aj 1 ' 1 A/jy + ... + N 1 H *Rki 2 

This equation looks like Eq.(56) but cannot be solved by the 
previously described iterations which involve refinement of 
weights since a, are centered at a fixed point xo(74), yo(Ti). 
However, a simple truncation of second-order terms A[’° Ap x 
and A^Apy and substitution of the measured A/i x for the true 
Ap x displacement yields 

An x = N~ l Y HTi)[xo(.T 2 ) - HTi )] - N~ l Y %2 (61) 


and a similar W y component of the vector IT. Note that W x is 
not zero as could be expected from (58) since positions and a, 
in Eq.(59) refer to different time moments. 

Definition (59) for W x looks like Eq.(19) except for the 
use of dj instead of optimal weights a, at which M x = 0 
and M v = 0. Therefore a function Q(q) corresponding to IT 
contains a small additional q 2 term which emerges since the 
first coordinate moments M x - H d/(7’i)[xo(7’ 2 ) - x;(74)] and 
M y = H d;(7’i)[yo(74) - y,-(74)] formed with the biased a I (7’i) 
values are not zero. 


The error caused by neglecting the difference between A/j x and 
A/r v is small and has a variance equal to that of image motion 
in a double star system with .v 1 = A// v ; it can be safely ignored 
for small A//. For instance, for distances A// = 4 mas measured 
with a D - 10 m aperture, any k. v - 3 and 7’ = 10 min the 
bias is about 340" • lO -6 y/Ap/40" ~ 3.4 /(as. The estimate was 
found by scaling the value of A given in Table 6 for a double 
star, from ,V| = 40" to the length A/( with a dependency A 2 ~ S 2 
. Even for such a large displacement, the relative error is about 
10 ~ 2 . 
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The reduction technique considered in this Section thus en¬ 
sures extracting of proper motions to within A accuracy, assum¬ 
ing of course that there are no other sources of noise. 

8. Image motion integrated over the atmosphere 

8.1. The model of C 2 vertical profile 

Because a total variance A 2 of image motion is equal to a sum 
of A 2 (/z) additives generated by each turbulent layer, its value 
therefore is a function of the C 2 (/z) vertical profile. In Fig.10 we 
reproduce typical plots of averaged C 2 (/z) for the three Chilean 
sites: Cerro Tololo ( http://www.gemini.edu/sciops/instruments/ 
adaptiveOpticsfSeeing.html ), Cerro Paranal (average over all 
profiles in Fig.2 given by Loaurn et al. 2000) and San Pedro 
Martir (Avila & Vernin 1998; average for the 1.5 and 2.1 m 
telescopes). The model profile by Hufnagel (1970) is shown 
for comparison. Since the data for h > 20 km are often un¬ 
available in original papers, we extrapolated C 2 up to h = 25 
km using a scaled profile from Serro Paranal (measured to 25 
km height) matched to the measured data at h - 20 km for the 
other sites. 

The divergence between local averaged profiles in Fig.10 
is about a half of a decade. However, the temporal variations 
of C 2 (/z) local profiles are much stronger, which Loaurn et al. 
(2000) had demonstrated in Fig.2 for the Cerro Paranal site. 
For these reasons it seems impossible to suggest a universal 
model of C 2 (/z) that will adequately match the real shape of 
the turbulence profile for any atmospheric conditions even at a 
single place. 

To derive numerical estimates, we defined the model of 
C 2 (/z) as an average of San Pedro and Cerro Tololo data which 
represent high and low limits of C 2 at /z > 20 km. Of course, 
the adopted model of C 2 (/z) is somewhat arbitrary and, due to 
varying atmospheric conditions, may give a factor 3-5 incor¬ 
rect predictions for a sample value of A. The model neverthe¬ 
less gives quite reliable estimates of average A for the Chilean 
sites and thus sufficiently well serves for the purpose of this dis¬ 
cussion. The vertical profile of the wind velocity V(/z) (Table 
5) which represents mean conditions for the South Geminy 
Telescope (Cerro Pachon, Chile) was taken from Avila et al. 
( 2001 ). 

Table 5. Wind velocity model taken from Avila et al. (2001) 


h, km 

V, m/s 

< 7 

10 

7-12 

20 

12-17 

40 

17-19 

20 

> 19 

10 


8.2. Contribution from different altitudes 

For narrow fields p/z <K D/ 2, application of the current method 
of differential observations is very promising since A goes now 



h, [km] 

Fig. 10. Average C 2 (/z) profiles for Cerro Tololo (1), Cerro Paranal 
(2), San Pedro Martir (3), Hufnagel's model (4) and the present model 
(squares) 

as a power k/2 (if k < v + p) or (v + p)/2 (if k > v + p) of the 
small quantity p/z/Z). For a 100 m telescope, the narrow field 
condition holds at any /z < 30 km providing the effective size 
p does not exceed 1 - 2'. For this reason, any increase of k 
and v parameters, corresponding to movement down the Table 
1 diagonal, always results in a better suppression of turbulent 
effects, especially those generated at low altitudes. In the case 
of a 10 m telescope, the choice of p is critical for the validity 
of the narrow field condition as even at quite moderate p = 1' 
it turns to be violated already at about h - 15 km. For upper 
layers, a much flatter, less efficient power dependency with an 
index p/2 holds, signalling a turn to a wide field mode of dif¬ 
ferential measurements (Lazorenko 2002a). 



h, km 

Fig. 11. A contribution A (h) from each 1 km thick atmospheric layer 
as a function of h, k and v for a 100 m telescope, 10 min exposure and 
1' effective field radius 

The model of C 2 (/z) defined in the above subsection was 
used to compute the contribution A(/z) from each A/z = 1 km 
turbulent layer to the total value of A 2 . The plots in Fig.-s 11 
and 12 display the function A(/z) computed for a 100 and 10 m 
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Fig. 12. The same as in Fig. 11 for a 10 m telescope 


telescopes placed at Iiq = 2.5 km altitude, T - 10 min, p = 2/3 
and stellar reference groups ”a, b, d, g, f' of Table 2, with k 
ranging from 2 to 12. The apodization parameter was near op¬ 
timal v = k + 1. Compare the effect of high order symmetriza- 
tion to the simple k = 2 (Fig. 11); for instance, with k — 8 and 
altitudes h < 10 km the value of A (h) is 5 orders lower than that 
with k = 2. At h ~ 20 km the gain still is 3 decades if D — 100 
m. 

While a contribution from low-altitude layers is large for 
asymmetric groups ( k = 2), symmetric groups with k > 4 or¬ 
ders are less sensitive to the turbulence at h < 15 km. The in¬ 
tegrated value of A for high k depends largely on high altitude 
turbulence in spite of the fast decrease of C 2 with h. For k > 8 , 
the function A (h) increases until h — 25 km, at least for the 
C\ profile adopted. The data on C 2 n behaviour at h ~ 25-40 km 
therefore is of special interest for the correct prediction of the 
image motion variance. Probably, the estimates of A computed 
in this study for k > 8 are slightly underestimated due to absent 
data for h > 25 km. 

8.3. The integrated variance of image motion 

Table 6 , similar to Table 1, gives values of A integrated over 
the atmosphere (in zenith direction) for apertures D — 4, 10, 30 
and 100 m. Calculations have been performed with the C, ](h) 
and V(h) model described above, p = 2/3, T = 10 min and 
telescope altitude ho = 2.6 km. The k values ranging from 2 to 
12 refer to tutorial configurations of Table 2, except the exam¬ 
ple of k = 4 order at entry ”c.f.” that presents a non-discrete 
stellar field with the Q(q ) function (48). Since the effective an¬ 
gular size p = V of each group corresponds to S » 4 - 6 m at 
15-20 km height, the condition of very narrow field S D/2 
is met only for D = 100 m, and, partially, for a D = 30 m. 

The estimates given for a 100 m telescope confirm the ef¬ 
ficiency of the high k and v used. In comparison to normal ob¬ 
servations (k = 2, v = 3), the gain in A is about 5 orders of 
magnitude with extreme parameter values given in Table 6 . The 
improvement in A which occurs with an increase of k at any v is 
also typical. On the contrary, the increase of v at fixed k is use¬ 
ful only up to v ~ k. A similar dependency is valid for smaller 


Table 6. The variance A (pas ) integrated over the atmosphere for star 
configurations ”a-h” of Table 1 and continuous field ”c.f.” for some k 
and v parameters; T = 10 min; p = V 



\ k 


2 

4 

4 

4 

4 

6 

8 

12 





conf 

. name 





a 

b 

c 

f 

c.f. 

d 

g 

h 

V 




D = 

= 4 m 




3 

572 

261 

244 

254 

233 

167 

114 

44 

5 

634 

315 

293 

305 

276 

219 

159 

58 

7 

679 

356 

330 

344 

307 

263 

200 

83 

9 

710 

387 

359 

374 

332 

297 

235 

106 

V 




D = 

10 m 




3 

337 

89 

87 

88 

85 

40 

25 

10.6 

5 

382 

114 

110 

112 

109 

46 

21 

4.9 

7 

419 

139 

134 

137 

132 

63 

29 

4.6 

9 

448 

162 

155 

159 

151 

80 

41 

7.0 

V 




D = 

30 m 




3 

165 

19 

19 

19 

19 

7.4 

4.80 

2.03 

5 

188 

21 

22 

21 

22 

3.9 

1.27 

0.26 

7 

207 

27 

27 

27 

27 

4.8 

0.98 

0.08 

9 

224 

32 

33 

32 

33 

6.4 

1.30 

0.05 

V 




D = 

100 m 




3 

74 

3.5 

3.5 

3.5 

3.4 

1.21 

0.783 

0.3038 

5 

84 

2.9 

2.9 

2.9 

3.0 

0.21 

0.061 

0.0111 

7 

92 

3.6 

3.7 

3.6 

3.8 

0.20 

0.017 

0.0010 

9 

99 

4.5 

4.6 

4.5 

4.7 

0.27 

0.018 

0.0002 

N min 

1 

3 

3 

3 

oo 

6 

10 

21 


apertures, with a tendency of the optimal v value to decrease to 
v — k — 1 or even v = 3 (no apodization) for a 4 m telescope. 
Note the very fast increase of A at transitions to smaller aper¬ 
tures. Nevertheless, considering entry for a 10 m telescope, one 
can note a progressive improvement in A with implementation 
of high k. Even with no apodization, the atmospheric error can 
be reduced to 10 p as providing that at least 21 reference stars 
are available in a circle of V radius (see Eq.(31)). Even limited 
to 6 stars, which does not allow k larger than 6, one can expect 
still quite small (~ 40 /.(as) errors suitable for exoplanet search 
programs. 

Symmetric ”b, f”, strongly asymmetric ”c” and non¬ 
discrete ”c.f.” configurations of equal k — 4 order have been 
purposely included in Table 6 to show that the magnitude of 
A depends rather weakly on peculiar features and type of star 
distribution in the field, and is, in fact, a function of k and v 
providing that p is fixed. 

Table 7 represents estimates of A for a reference field given 
in Table 3, at 10 min exposure; the last line contains effective 
sizes p for each k. Note that a good quasi-symmetric distribu¬ 
tion of stars alleviates the difference in A between k — 2 and 
k = 4 orders, the noise for k - 2 is only slightly over that given 
for k = 4. 

Pravdo & Shaklan (1996) derived an estimated A = 150 
/(as/hr for the magnitude of atmospheric fluctuations at the 5 
m telescope and Table 3 star field. This is a value obtained by 
considering each field star, in turn, as a target, and by averaging 
individual estimates of A that fluctuated at least a factor of 2- 
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Table 7. The variance A (//as) for the star field of Table 3 at 10 min 
exposure 



1 k 


2 

4 

6 

8 

V 


D = 

: 10 m 


3 

60 

53 

14 

9.0 

5 

75 

67 

13 

5.8 

7 

91 

82 

18 

7.8 

9 

105 

96 

23 

11.1 

V 


D = 

100 m 


3 

5.1 

2.1 

0.420 

0.268 

5 

5.5 

1.7 

0.050 

0.016 

7 

6.2 

2.2 

0.042 

0.004 

9 

6.8 

2.6 

0.056 

0.004 

P 

3.7" 

45.7" 

35.2" 

42.2" 


3. Following the target position changes, the effective frame 
size (approximately equal to the frame radius with the target at 
its center) varied from 3.7" to 90" for the outer stars. For this 
range of p, our model predicts variations from 55 to 300 //as/hr 
which well matches the observed value of 150 //as/hr. 

9. Astrometric performance of very large 
ground-based telescopes 

Realization of 1-10 //as accuracy requires a good elimination 
of various noise sources related to optical aberrations, pixeliza- 
tion effects (especially for small images produced by adaptive 
telescopes), photon noise in star images, differential chromatic 
refraction (DCR) etc. As it was noted by Louarn et al. (2000), 
in particular, the problems caused by a DCR that stretches the 
star images into colored strips are very intricate. The ampli¬ 
tude of the DCR effect depends on zenith distance, air temper¬ 
ature and pressure, spectral band and star colors. Using Allen’s 
(1973) tables, one can find that two rays with wavelengths of 
500 and 600 nm coming from a star at a zenith distance of 30° 
are imaged with a separation of about 180 mas along a verti¬ 
cal direction. The noise induced by this effect in the differential 
position of stars (Pravdo & Shaklan 1996) amounts to about 60 
/./as in a 1.5' field for the 5-m Palomar telescope. Fortunately, 
in proper motion studies the DCR effect is residual and essen¬ 
tially weakened since star motions are found from residuals of 
differential star positions in the two epochs. The modelling of 
DCR based on use of atmospheric bulk parameters is there¬ 
fore very promising. We have found that a proper control of 
atmospheric air parameters (air temperature to 0.2°, pressure 
to 0.2 mb) allows one to apply corrections which reduce DCR 
noise to 8 //as in the relative displacement of A and M stars 
in the field of 1.5'. Once effective wavelengthes of stars are 
known to 0.4 nm, the noise decreases to 0.8 //as. It should be 
noted that for high quality adaptive optics producing images 
with FWHM<s 100 mas, the DCR corruption of images is so 
strong that it makes them entirely unsuitable for measurements. 
It is necessary therefore to use some special optics for compen¬ 
sation of atmospheric chromatism, otherwise the filter width 
should be strongly narrowed. 


Assuming that solution of this and other problems will be 
eventually found by progresses in technology, we restrict the 
error budget with two components: the atmospheric image mo¬ 
tion and photon noise with variances A 2 and cr 2 ph respectively. 
The contribution from both effects was evaluated as a function 
of the angular field size R for sky star densities near the galac¬ 
tic plane and at the pole. All particular cases of aperture, image 
parameters, exposure, field size etc. of course, not can be con¬ 
sidered; therefore we restricted analysis only to the case of a 
future extremely large 100 m telescope and modern 10 m class 
telescopes. 

Estimates of cr p h were found with use of Eq.(47). The value 
of FWHM in this expression strongly depends on the perfor¬ 
mance of adaptive optics, the telescope aperture and may vary 
from 0.0015" (diffraction limit of a 100 m telescope) to 0.4" 
(atmospheric uncorrected seeing). The next estimates for a 100 
m telescope assume FWHM=0.1" achievable with low-order 
adaptive optics, and for a 10 m telescope a quite conservative 
FWHM=0.4" was adopted. We assumed that observations are 
obtained in zenith, in R band, CCD quantum efficiency 0.85, 
transmission of optics 0.8 and of atmosphere 0.9. Then a star 
of V=15 mag and of average spectral type, being observed 
with a 100 m telescope, will provide 0.93 • 10 s detected elec¬ 
trons/sec (Allen 1973). A total light n t of the star field was 
estimated based on the Galaxy model by Bahcall & Soneira 
(1980). Stars fainter than V=23 mag were not considered as 
they give low light signal. To obtain more robust results, the ex¬ 
pected star number in each 1 mag bin was rounded (truncated) 
to the smaller integer. This procedure trimmed the bright end 
of stellar magnitudes due to which very narrow star fields were 
formed largely by the faintest stars. 

Fig.13 represents data for a 100 m telescope, T - 10 min, 
galactic coordinates b = 20 °, l = 0 ° and k in the range from 
2 to 12. For k < 8 , the apodization parameter was set to be 
v = k + 1, for higher k its value was limited by v - 9 so as not to 
worsten the light transmission. The dashed lines represent cr p h 
computed with Eq.(47); plots start from the smallest field size 
which ensures N mln star number necessary to realize k order 
symmetry. 

The estimates of A were computed by integration of Eq.(12) 
for the turbulence model described in Section 8 . It was per¬ 
formed two ways: 

1) by processing random stellar fields simulated with a Galaxy 
model (Bahcall & Soneira 1980). It involved calculation of a, 
defined by conditions (28), (32) and of Q(q) filter function (24) 
for each field. These direct point estimates A nl( | of atmospheric 
error A computed for R < 0.8" are shown by dots for odd k/2 
and by crosses for k /2 even. 

2 ) by computing the atmospheric error for ’’typical” stellar 
fields whose filter Q(q ) mathematical expectation is given by 
Eq.(52) for odd k/2 (parameters y and p taken from Table 4) 
and (48) for even fc/2. Computed estimates of A exp are shown 
by solid lines. 

The difference between A ln( j and A exp estimates is caused 
by different values of sampled effective sizes S rn( i and of S exp 
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Fig. 13. Atmospheric and photon noise as a function of field radius 
R for a 100 m telescope, 10 min exposure, near the Galactic equato¬ 
rial plane and a set of k and v parameters. Atmospheric noise: dots 
- for random star fields (A rn j) at k/2 odd; crosses - the same at k/2 
even; solid lines - expected mean noise A exp for "typical” stellar fields; 
circles - transformation (62) from A md to A exp (see text); triangles - 
reversed transformation. Dashed lines - photon noise cr ph at 0.1" see¬ 
ing 

(the expected effective size). The inter-relation of these values 
is given by the expression 

Aexp = A m d(S exp /s md) k/2 (62) 

that follows from the power law (38) valid for k < v + p. Open 
circles in Fig. 13, the results of transformation (62) from A rn d 
to A exp , are shown to be placed perfectly along solid lines that 
represent A exp . 

The expression opposite to Eq.(62) can be used for indirect 
computation of A rn( j proceeding from an effective frame size 
S m( | and A exp for fields containing N ~ 10 2 -10 3 stars. Direct 
numeric integration in this case is too time-expensive since the 
number of terms in Eq.(24) increases as N 2 or /( 4 . With this 
approach, atmospheric noise is easily estimated at any large R. 
Estimates of A rn d computed based on S m d and A exp are shown 
in Fig. 13 by triangles for k = 8, v = 9 and R varying from 3 
to 8' when star number in the field mounts from 1000 to 3000. 
For R < 0.8', approximate estimates (triangles) exactly match 
those directly computed A rn d (crosses). 

The positive factor 2 offset of A exp over point estimates A rn d 
seen for k = 4,8,12 and very narrow R < 0.6' fields with 


N < 20, originates from use of a non-descrete field model for 
computation of A exp , or the assumption S — R (50). At very low 
N , however, this gives rather an upper limit of 5 but its mathe¬ 
matical expectation since stars do not entirely cover periphery 
of the field. The difference is small but becomes apparent after 
being amplified due to the power dependency A ~ S k ^ 2 . 

A short comment should be made concerning plots for 
k = 10 and k = 12 with the noticeable scattering of A exp (cir¬ 
cles) computed with Eq.(62). For these plots, observations are 
carried out under condition k > v + p, which, according to 
Eq.(36), means that A 2 is proportional to the power v + p of 
modified frame size S which is not equal to S. Though Eq.(62) 
is not valid here, it still gives rather good results. 



Fig. 14. The same as in Fig.13 (D = 100 m), for the Galactic pole; 
data for random star fields are not shown 


A plot with error estimates for the Galactic pole is given in 
Fig. 14. For a 10 m telescope, results are represented in Fig.- 
s 15, 16; no apodization is applied since it does not offer an 
improvement (Table 6 ). 

A peculiar feature of Fig.-s 13-16 is the critical point Rq 
at which plots of and A drawn for a certain k and v inter¬ 
sect. For R < Rq, where filtration of atmospheric noise is very 
efficient, photon noise dominates over A (photon-limited ob¬ 
servations). For R > Rq, on the contrary, A > cr p h and observa¬ 
tions are atmosphere limited. At low sky star density, reference 
frames containing at least N m j„ stars (brighter than V=23), nec¬ 
essary to form k order symmetry fields, can be formed only at 
relatively wide R where A becomes large. The typical situation 
described is shown in Fig. 16 (polar regions) where plots for 
cr p h and A expected for a 10 m telescope do not intersect at any 
R, signifying that observations are atmospherically limited. 

Consider now a total error cr t which we define simply as 
being equal to the largest error component cr p i, or A. At any R, 
a value of cr, can be minimized by proper selection of k and 
v. The plot of optimal cr, error as a function of R was formed 
taking the best combination of segments of cr ph and A curves in 
Fig.-s 13-15. Resulting segmented curves drawn for the cases 
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field radius, arcmin 


Fig. 15. The same as in Fig. 13 (Galactic equator), for a 10 m tele¬ 
scope and 0.4" seeing; data for random star fields are not shown. No 
apodization (v = 3) 



1 10 
field radius, arcmin 

Fig. 16. The same as in Fig. 15 (a 10 m telescope) for a Galactic pole 


discussed are shown in Fig. 17. Each curve starts at small R 
where only low order k = 2 or k — 4 symmetry is realized 
with 1-3 reference stars; the right ends of curves correspond to 
k — 12 and fields containing N ~ 100-1000 stars (except for 
the case of a D - 10 m telescope operating in the polar region 
where N ~ Nu m ). Plots show that cr, is a decreasing function 
of R , so use of high k > 12 orders results in a progressive slow 
decrease of cr,, which, however, involves the use of a wider 
field. From Fig.17 data we may notice that the near optimal 
field size is 0.4-1' at high star density and about 2' at polar 
regions. 

With above assumptions on FWHM, and at optimal field 
size, the expected error of ground-based observations for a 10 
m telescope (no apodization) is, depending on sky star density. 



field radius, arcmin 

Fig. 17. Total optimized atmospheric and photon noise error cr, as a 
function of field radius R for 10 and 100 m telescopes, a 10 min expo¬ 
sure and parameters of observations same as in Fig.-s 13-16; solid 
lines — for Galactic equator (b = 20°, l = 0°); dashed - for the 
Galactic pole. Along each line, k value increases with R from 2 to 
12 

10 to 60 //as per 10 min exposure. For a 100 m telescope this 
estimate is 0.2 to 2 //as. 

10. Conclusion 

As commented by Findegren (1980), dependency A ~ 
(ph/D)D l ^T~ 1 ^ 2 valid for conventional differential astrometric 
observations in narrow fields (in wide fields A ~ (ph) l ^ 3 T^ 1 ^ 2 ) 
is not a fundamental limitation to observations from the ground 
as it can be improved by applying a new method of measure¬ 
ments. In this paper we introduce a method that essentially im¬ 
pairs restristions on the precision of ground-based astrometry 
caused by atmospheric turbulence; the power dependency of 
A on ph, D and T now becomes much stronger (53). Efficient 
filtration of atmospheric wave-front distortions is achieved pri¬ 
marily due to application of reference fields of enhanced vir¬ 
tual symmetry. The method takes advantage of using very large 
telescopes. Thus, the total error of observations (atmospheric 
and photon noise) for a 10 m telescope and non-corrected im¬ 
ages is expected to be about 10-20 //as per 10 min exposure, 
providing other sources of errors are small. In some cases, at 
least, for the study of extra-solar planets, such an accuracy is 
acceptable. Of course, very high precision ground-based obser¬ 
vations can be performed only in narrow fields; global data are 
to be obtained from space. We hope the method described is 
applicable in practice and that the results of this study will be 
of interest. 
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